Resonating valence bond wavefunctions and classical interacting dimer models 
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We relate properties of nearest-neighbour resonating valence bond (nnRVB) wavefunctions for 
SU{g) spin systems on two dimensional bipartite lattices to those of fully-packed classical dimer 
models with potential energy V on the same lattice. We define a cluster expansion of V in terms 
of n-body potentials Vn, which are recursively determined from the nnRVB wavefunction on finite 
subgraphs of the original lattice. The magnitude of the n-body interaction V„ (n > 1) is of order 
<D{g~^"^^^) for small g~^, while Vi reduces to a constant due to the fully-packed nature of the model. 
At leading non-trivial order on the square lattice, the interacting dimer model only has two-body 
interactions V2 (5) that favour two parallel dimers on elementary plaquettes. Setting g = 2 and using 
the results of earlier work on this interacting dimer model, we find that the long-distance behaviour 
of the bond-energy correlation function is dominated by an oscillatory term that decays as 
with a ~ 1.22 for SU{2) spins. This result is in remarkable quantitative agreement with earlier 
direct numerical studies of the corresponding wavefunction, which give a « 1.20. 

PACS numbers: 75.10.Jm 



Spin liquid states of low dimensional insulating mag- 
nets, in which the constituent spins fail to develop any 
kind of long range order down to T = in spite of strong 
magnetic exchange interactions, are an interesting con- 
sequence of competition between the different magnetic 
interactions in the system. Explicit variational wavefunc- 
tions that encode such behaviours have played a very 
important role in our understanding of such states of 
matter Perhaps the best known examples of such 
wavefunctions are the resonating valence bond (RVB) 
wavefunctions introduced by Anderson and collabora- 
tors 

The simplest of these is the nearest neighbour RVB 
(nnRVB) wavefunction |'4'(g)) for SU{g) spins on a two 
dimensional bipartite lattice. R is written as a uniform 
amplitude superposition of all possible product states in 
which each A-sublattice spin forms a SU{g) singlet state 
with one of its _B-sublattice neighbours Q . Although the 
g = 2 wavefunction has been studied on the square lattice 
for over twenty years now we owe a detailed under- 
standiiig of its properties to much more recent work 043 '■ 
In Ref [ll,^], spin and bond-energy correlations were mea- 
sured in the square lattice case using Monte-Carlo meth- 
ods to establish that this state has an exponentially de- 
caying spin correlation function, \Cs(f)\ ~ exp( — jf^/^), 
but power-law bond-energy correlation functions at large 
1^1: \CEir)\ ~ with a « 1.20. While such a 

short-ranged spin correlation function is a characteris- 
tic property of spin-liquids, the power law form of the 
bond-energy correlation functions strongly suggests that 
the nnRVB wavefunction for SU (2) spins on the square 
lattice actually describes a critical state on the verge of a 
transition to an ordered phase in which the bond-energies 
order. 

Here, we develop a precise non-perturbative mapping 
that connects properties of \^{g)) on a two-dimensional 
bipartite lattice to those of a classical fully-packed dimer 



model on the same lattice, which has a non-trivial in- 
teraction potential V for the dimers in addition to the 
usual non-overlapping constraint. We define a cluster 
expansion of V in terms of n-body potentials Vn , which 
are recursively determined from the nnRVB wavefunction 
on finite subgraphs of the original lattice. Vi reduces to 
a constant due to the fully-packed nature of the dimer 
model, while the n-body interaction I4, (n > 1) is of or- 
der 0((7~("~^^) for small g~^, and thus decreases with 
n. The rate of decrease is controlled by the smallness 
of g~^, which also controls, via this mapping, the expo- 
nential decay of spin correlation functions in 1^(5)). To 
leading non-trivial order on the square lattice, the inter- 
acting dimer model only has two-body interactions V2 (g) 
that favour two parallel dimers on elementary plaquettes. 
Setting .9 = 2 and using the results of earlier work 
on this interacting dimer model, we find that the long- 
distance behaviour of the bond-energy correlation func- 
tion in |'4'(,9 = 2)) is dominated by an oscillatory term 
that decays as with 



1.22 



(1) 



This result is in remarkable agreement with recent 
studies of the SU{2) nnRVB wavefunction 0, @], and 
provides a quantitative resolution of the surprising co- 
existence of short-ranged spin correlations and power- 
law bond-energy correlations in \^{g = 2)). Our non- 
perturbative mapping to an interacting fully-packed clas- 
sical dimer model explains the success of the phenomeno- 
logical approach of Tang et. al. Q , who were motivated 
by the connection between fully-packed dimer configura- 
tions and ground states of certain [ll| large-A^ SU{N) 
square lattice antiferromagnets to develop a phenomeno- 
logical height model description of the critical aspects of 
\'^{g = 2)), similar to that used for the long-waveleiigth 
properties of classical fully-packed dimer models 12, ll3| . 
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We begin our analysis by writing 



= 5Z 125)9 '^here 
ma - n \^^^9))e ■ 



(2) 
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Here, the summation is over all fully-packed dimer con- 
figurations D, the product is over all edges e covered by 
a dimer in 2), and |$o(.9))e is the SU{g) singlet state of 
the two spins connected by edge e. The norm 



{^{g)\^{g))= Y.i'D'm, 



(3) 



'D,D' 



can be readily expressed as the partition function Z(g) 
of a classical fully-packed loop model in which each site 
belongs to exactly one non-intersecting loop, which can 
either be a doubled edge (consisting of a single edge tra- 
versed in both directions) or a non-trivial loop (consisting 
of four or more distinct edges which are each traversed 
once). For the case of SU{g) spins, one obtains 



(4) 



with wioopig,^) = (g)"<'(^)(2g)"'(^', where nd{L) is the 
number of doubled edges and ni{L) the number of non- 
trivial loops in the loop configuration L (Fig[T]). Expec- 
tation values and correlators of various physical quan- 
tities can also be represented in terms of estimators 
that measure probabilities for various loop-gas configu- 
rations. In the SU (2) case, the spin correlation function 
Cs(r) = ('^'(O) • 'S!{'r)) is proportional to the probabil- 
ity that points and r both lie on the same loop. On 
the other hand, the connected bond-energy correlation 
function C£^(f) = (B^(O)B^(r)) - (5^(0)) (B^(f)) can 
be related to the probabilities for the four points 0, fi, 
r, and r -\- fi. to lie on the same loop or at most on two 
different loops Here, B^(f) = (5(r) • S{r-\- fi)) with 
(1 ~ x,y representing an elementary lattice translation 
in the x or y direction. Since Csir) is proportional to 
the probability of points and r being on the same loop, 
the short range nature of Cs implies that the g = 2 loop 
model is in a "gapped" phase with predominantly short 
loops. Indeed, since Zioop defines a conformally invariant 
loop model with a power law distribution of loop sizes 
for g = 1 [l3|j we expect that this distribution becomes 
exponential for g > I. The most natural scenario is then 
that the entire g > 1 short-loop phase of Zioop is con- 
trolled (from a renormalization group standpoint) by the 
g = 00 fixed point. 

Our approach to the existence of power law bond- 
energy correlation functions in this short loop phase 
starts with the elementary observation that a loop gas 
with extremely short loops can still have long-ranged cor- 
relations in the position and orientation of loops. More 
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FIG. 1: (2) |I')g can be represented in terms of the loop con- 
figuration £ obtained from the overlap loops in the overlap 
diagram generated from dimers in T) and D . Each overlap 
loop contributes a factor of g to {D \D)g due to the overlap 
of the SU{g) singlet states that make up |D )g and \'D)g. By 
convention, all overlap diagrams related by independent inter- 
change of black and shaded dimers in each overlap loop are 
identified with the same loop configuration L, and therefore 



precisely, at g = 00, all loops are trivial in that they have 
length 2, and correspond to doubled-edges. The weight 
of any such loop configuration made up entirely of dou- 
bled edges is always (7^=/^ (where Ns is the number of 
lattice sites). We now represent all these doubled-edges 
by dimers to map such a loop configuration to a fully- 
packed dimer configuration T) on the same lattice. Thus, 
loop configurations that survive the g ^ 00 limit all have 
equal weight and define the partition function of a fully- 
packed dimer model in which each dimer contributes a 
factor of g to the weight of D: 
independent of D. 



For g < oo, 2,/oop also gets contributions from more 
general configurations £ consisting of both doubled-edges 
and non-trivial loops (with 4 or more edges). Each non- 
trivial loop in such a finite-g configuration L can be re- 
placed in exactly two ways by a sequence of doubled- 
edges on alternating edges of this non-trivial loop. Thus, 
a general loop configuration Z with ni non-trivial loops 
and Ud doubled edges maps to 2"'^''"^ different loop con- 
figurations made up purely of doubled edges, which we 
represent by dimers. Each finite-g configuration L of 
the original loop model thus maps to 2"'^''"^ dimer con- 
figurations D„ (a = 1,2...2"'('^)). Next, we distribute 
wioopig, the original weight of £, equally among these 
Da. As a result, each of these 2"''^''') different configura- 
tions "Da acquire an additional weight w{g, L)/2"'^''''K 

This maps the original loop model to a dimer model 
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with weights 



Wdi'. 



{^{g)\V),, (5) 



where L\T> denotes all loop configurations L compati- 
ble with the fully- packed dimer cover D, i.e. obtain- 
able by superposing some fully-packed dimer cover T> 
on the given V, and we have used Eqn ( [3]) and ^ 
to obtain the second equality. The original loop par- 
tition function 2,;oop is then equal to the partition sum 
over all fully-packed dimer configurations 2) with weights 
(5,2)): 



loop 



dime 



Wdi: 



(6) 



As Wdimerig,'^) depends on the structure of 2), we 
have provided a precise non-perturbative mapping of the 
original loop model to an interacting, classical, fully- 
packed dimer model with potential energy V{g, D) given 



(7) 



We now define a decomposition of the potential energy 
V{g, 2)) of a fully packed dimer configuration D into a 
sum of n-body potential energies Ki (D„) of subconfigu- 
rations D„ consisting of n distinct dimers from D: 



(8) 



n D„G'D 



The Vn are determined recursively from computation 
of the weight w^ij^^rid^'^") '^^ 2?ti in the interacting 
dimer model on the finite subgraph Qni'^n) of the square 
lattice. Here, this weight is calculated via Eqn ([5|) from 
the loop model defined on 9ri(25n), and the subgraph 
9n(2)ri) consists of the 2n vertices covered by dimers of 
the subconfiguration Dn, along with all allowed edges 
between these vertices. 

In the first step of this recursive construction, we 
consider any particular Di and determine the weight 
w dimer id ^'^i)- The Original loop model on Si('Di) has 
only one valid configuration, which is a doubled edge on 
the only edge of Si(Di). Using Eqn ([5|), we therefore 
have w dimer (9 t'^'^) ~ 9- "^^^ logarithm of this weight 
defines the one-body potential 



(_) = V^ (I) 



log(5) ■ 



(9) 



With this in hand, Vn for arbitrary n can be obtained 
recursively from the equation 



los 



dimer 



n-1 



m=l D„eD„ 

(10) 



where Dm S Dn denotes all m-dimer subconfigurations 
Dm of Dn, and 9m{Dm) the corresponding subgraphs of 
S„(D„). 

For instance, to obtain V2, we consider any particu- 
lar D2 and determine WdimerC^^) follows: If 32(2)2) 
does not form a plaquette of the square lattice, there is 
only one valid configuration of the loop model on 32(2)2) 
(consisting of doubled edges that take the place of the 
dimers in 2)2). From Eqn (P, this gives w^^^^^{D2) = 5^. 
On the other hand, when 32(2)2) does form a plaque- 
tte of the square lattice, w^imeri'^'^) S^ts contributions 
from two of the three valid configurations of the loop 
model on 32(2)2)- One of these consists of two doubled 
edges that take the place of the dimers in D2, while the 
other is a non-trivial loop on the boundary of the plaque- 
tte. Eqn (O then gives w^-^^^{D2) = .g^ -f g. Knowing 
^dimer(2^2), V2 Can bc obtained from the recursion rela- 
tion Eqn pI7|) with n = 2. Clearly, V2 is non-zero only if 
the two dimers live on the same plaquette of the square 
lattice, and in this non-trivial case we obtain 



1/2 (□)= 1^2(0) = -log(l + .g-i) 



(111 



Similarly, V3 is seen to be zero unless the three dimers 
live on a pair of adjacent plaquettes, and in this non- 
trivial case we obtain 

V3pj)=-log[{l + g-'+g-'){l + g-Y'] , 
V, (m) = - log [(1 + 2g-'){l + g-Y^] , (12) 

and the symmetry-related counterparts of Eqn (jl2p ob- 
tained by using lattice reflection and rotation symme- 
tries. 

It is clear that each Vn for n > 1 is of order (3{g~'^^~^'>) 
when the n dimers live on a contiguous set of plaquettes 
of the square lattice (sharing edges with each other), and 
zero otherwise. Our procedure thus expands V in the 
size of clusters, and is controlled by the smallness of g"^. 
It may be viewed as a classical version of the "Schricffer- 
Wolff " canonical transformation approach familiar in the 
context of strongly interacting electronic systems; in the 
usual [l3l Schrieffer- Wolff transformation, the effects of 
higher energy states in a larger Hilbert space are encoded 
in modifications to the effective Hamiltonian that acts in 
the subspace of low-energy states, while in our classi- 
cal version of this approach, the effects of lower weight 
configurations of the full loop model are encoded in the 
effective weights of a restricted subclass of high-weight 
dimer model configurations D. 

Expectation values of any operator O can be calculated 
using the modified estimator (2)) constructed from the 
original loop gas estimator [IJ, ll6| J'q(£) to correctly 
encode the effects of the low-weight configurations that 
are not contained in the dimer model: 



5'o(2)) 



Wdv. 



,D) ^ 



UJloop 

2»i('C) 



(13) 
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Consider for instance the spin correlation function. To 
zeroth order in g~^, the corresponding modified estima- 
tor is non-zero only if and f are connected by a dimer in 
D, while the 0{g~^) correction gives a non-zero result if 
and r belong to the same flippable plaquette of 2) with- 
out being connected by a dimer, and similarly for higher 
order corrections. Thus the spin correlation function is 
expected to decay exponentially as a function of spatial 
separation. In the case of the bond-energy correlation 
function CE^j,{r) defined earlier, it is readily seen that 
the leading large \r\ behaviour is dominated by the first 
0(g°) term in the modified estimator, which is obtained 
by replacing each B^j^ (r) by the dimer occupation number 
n^(r) of the corresponding edge in the dimer model: 

CE^r^ ~ (n,,(0)n^(f))y - (n^(0))y (n^lr'))^ , (14) 

where {■■■)¥ denotes averages computed in the dimer 
model with energy V . 

As is well-known, such fully-packed interacting dimer 
models admit a microscopic height representation liol . 
121, ll3|, which upon coarse-graining leads to a coarse- 
grained height action that takes the form 



S ^ Tip I d^r{Vhf + 

p=4,8,12. 



Vp 



j (Prcos(2Trph) 



(15) 



where the ellipses denote higher gradient terms and 
highe r p owers of gradients consistent with symme- 
tries 
of g. 



12| , and the bare coefficients p and yp are functions 



The renormalization group theory for this height 



action is standard [12|, |19|: In the present variables, it 
tells us that there is a line of fixed points p* = K,y* ~ 
with < K < 4. As long as the bare values of p, yp 
and the coefficients of the omitted higher derivative and 
nonlinear terms are not too large, the system flows to an 
attractive fixed point K(g) on this fixed line. 

As mentioned earlier, the g — > oo limit maps to a non- 
interacting dimer model since all interactions {Vn with 
n > 1) vanish and Vi simply represents the fugacity 
g of each dimer. Therefore expect K{g — > oo) = 1/2 
since this is the known value of the stiffness for a non- 
interacting dimer model on the square lattice [l2j. As g 
is reduced from g = cx), the leading effect is an interac- 
tion V2 that favours flippable plaquettes. The interact- 
ing dimer model with only V2 present has been studied 
in detail in Ref , which established that the renor- 

malized stiffness k increases monotonically with the mag- 
nitude of V2 until it reaches k = 4, at which point the 
system undergoes a Kosterlitz-Thouless like transition to 
a columnar ordered state. 

To access the physics of the nnRVB wavefunction of 
SU{2) spins, we therefore identify the magnitude of V2 
with the inverse temperature (3 of Ref M, and set g = 2 
to obtain /3 = log(l + g'^) « 0.405 This places 



us deep in the high-temperature phase well above the 
transition to columnar order, and from Fig. 31 of Ref 0], 
we obtain the estimate K{g = 2) « 0.82. Calculating 
the required dimer correlation function from the fixed 
point height action with this value of k using the standard 
correspondence between dimer occupation numbers and 
the height-field [ill, [l^, we obtain the leading, large \r\ 
form of the bond-energy correlation functions: 



|f|l/K(9) 



|f|i/«(ff) 



(16) 



which gives, upon setting g = 2, the result advertised ear- 
lier. In addition, both these correlation functions have a 
subdominant piece which goes as (—1)^+^/1?^^ indepen- 
dent of n{g). 

We conclude by noting that our approach generalizes 
in an obvious way to include singlets between A and B- 
sublatticc sites that are further apart on the square lat- 
tice, as well as to the SU{g) nnRVB wavefunction on the 
honeycomb lattice. Again, we expect the corresponding 
loop models to be in a short-loop phase for all g > 1, 
and this allows us to access the physics of these wave- 
functions for 3 > 2 by mapping to an interacting dimer 
model; on the honeycomb lattice, the leading interaction 
terms will favour flippable hexagons, while the presence 
of longer range singlets will introduce additional interac- 
tions on the square lattice. Possible generalizations to 
various three-dimensional bipartite lattices are more in- 
triguing, since it is not obvious that the loop model will 
be in a short-loop phase for all g > 1 on any of these 
lattices. 
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